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Sunyaev-Zel'dovich observation of the Bullet-like cluster Abell 2146 
with the Arcminute Microkelvin Imager * 
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ABSTRACT 

We present 1 3.9 -18.2-GHz observations of the Sunyaev-Zel'dovich (SZ) effect towards Abell 
2146 using the Arcminute Microkelvin Imager (AMI). The cluster is detected with a peak 
signal-to-noise ratio of 13<x in the radio source s ubtracted map from 9 hours of data. Compar- 
ison of the SZ image with the X-ray image from lRussell et a i] (120101) suggests that both have 
extended regions which lie approximately perpendicular to one another, with their emission 
peaks significantly displaced. These features indicate non-uniformities in the distributions 
of the gas temperature and pressure, and suggest complex dynamics indicative of a cluster 
merger. We use a fast, Bayesian cluster analysis to explore the high-dimensional parameter 
space of the cluster-plus-sources model to obtain robust cluster parameter estimates in the 
presence of radio point sources, receiver noise and primordial CMB anisotropy; despite the 
substantial radio emission from the direction of Abell 2146, the probability of SZ + CMB pri- 
mordial structure + radio sources + receiver noise to CMB + radio sources + receiver noise 
is 3 x 10 6 : 1. We compare the results from three different cluster models. Our preferred 
model exploits the observation that the gas fractions do not appear to vary greatly between 
clusters. Given the relative masses of the two merging systems in Abell 2146, the mean gas 
temperature can be deduced from the virial theorem (assuming all of the kinetic energy is 
in the form of internal gas energy) without being affected significantly by the merger event, 
provided the primary cluster was virialized before the merger. In this model we fit a simple 
spherical isothermal /J-model to our data, despite the inadequacy of this model for a merging 
system like Abell 2146, and assume the cluster follows the mass-temperature relation of a 
virialized, singular, isothermal sphere. We note that this model avoids inferring large-scale 
cluster parameters internal to r2oo under the widely used assumption of hydrostatic equilib- 
rium. We find that at r^oo the average total mass Mj = (4.1 ± 0.5) x 1O 14 /i _1 M and the mean 
gas temperature T = 4.5 + 0.5 keV. 



1 INTRODUCTION 



Galaxy clusters are the largest collapsed structures known to exist 
in the Universe. The masses of rich clusters can reach as 10 I5 /i _1 M Q 
and the more distant ones, from around z > 0.2, subtend several ar- 
cminutes on the sky due to the slow variation of the angular di- 
ameter distance with redshift. As a result, clusters are powerful 



tracers of structure formation and evolution on scales of the or- 
der of a few megaparsecs. According to the standard A Cold Dark 
Matter (ACDM) model, galaxy clusters form via hierarchical in- 
teractions of smaller subsystems. During merger, these subclusters 
collide at relative velocities of thousands of kms~' and can release 
gravitational binding energies of up to ~ 10 57 J, which can lead 
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to shocks in the intracluster medium (IM). These conditions make 
cluster mergers ideal places to study the dynamics of matter under 
extreme conditions. The three assumed main components compris- 
ing the cluster, namely galaxies, hot ionized gas and dark matter, 
exhibit very different behaviours during subcluster mergers. The 
hot intergalactic gas is heated and compressed by the hydrodynam- 
ical shocks produced during the passage of the subcluster through 
the core of the prima ry, whereas the dark ma tter and galaxies are 
collisionless (see e.g. lMarkevitch et al.l ( 120070 ). As a result, the gas 
is slowed down by ram pressure and is displaced from the dark 
matter and the galaxies. Later, when the subcluster reaches regions 
of lower gas density in the primary cluster, the ram pressure drops 
sharply. Without as much ram pressure, the gas pressure and sub- 
cluster gravity cause some of the subcluster gas, which had been 
lagging behind the subcluster's dark matter centre, to 'slingshot' 
past it. This gas is then left unbound from t he subcluster and free 
to expand adiabatically dHallman et al.l l2004). 

Abell 2146 is a cluster at z = 0.23 consisting of two merging 
subclusters. The smaller subcluster passed through the centre of the 
larger subcluster some 0.1 - 0.3 Gyrs ago producing shoc k fronts 
which have been detected by Chandra JRussell et alj|2010h . These 
shock fronts are unusual features which only show at a specific 
stage in the cluster merger, before the shock reaches the outer, low- 
surface-brightness regions, and at angles on the sky plane which 
usually prevents the projection from hiding the density edge. There- 
fore, it is not surprising that shock fronts with Mach numbers sig- 
nificantly greater tha n one have only been detected in two other 
clusters: 1 E0657-56 dMarkevitche t al. 2002) -the "Bullet cluster"- 
and A520 jMarkevitchll2006l) . Unlike A520 the Bullet cluster and 
Abell 2146 appear to be at an early stage of the merger event, where 
the cluster dynamics are simpler and the separation of the hot gas 
and the dark matter components is clearer. 

The thermal Sunyaev-Zeldovich (SZ) effect provides an inde- 
pendent way of exploring the physics of the intracluster gas and 
examining typical cluster parameters such as core radius and gas 
mass. When Cosmic Microwave Background (CMB) photons tra- 
verse a rich galaxy cluster some will be inverse-Compton scattered 
by t he random thermal motion of the electrons in the intracluster 
gas dSunvaev &Zel'dovichlll970l : lBirkinshaw|[l 999). Unlike X-ray 
surface brightness, SZ surface brightness is independent of redshift 
and is therefore well suited for the study of galaxy clusters at any 
redshift. It is also less sensitive than X-ray measurements to small- 
scale clumping and the complex dynamics associated with the clus- 
ter core. 

In this paper we present 16-GHz SZ effect images of Abell 
2146 using AMI. In Section[2]we discuss the telescope, while de- 
tails of the observations and the reduction pipeline are given in Sec- 
tion[3] In Section|4]Bayesian inference is introduced. Section[5]de- 
scribes our analysis methodology while in Section [6] and Section 
[7] we present the results and discuss their significance. We present 
our conclusions in Section[8] 

Throughout the paper we assume a concordance ACDM cos- 
mology with a m>0 = 0.3, n A-0 = 0.7, Q. k = 0, Q. b = 0.041, w = -1, 
w a = 0, erg = 0.8 and H = lOOkms^'Mpc -1 . Relevant param- 
eters are given in terms of the dimensionless Hubble parameter 
h = H /100 kms~' Mpc~', except where otherwise stated. We also 
refer to h x = H /X kms~' Mpc~'. All coordinates are at epoch 
J2000. 



Table 1. AMI technical summary 





SA 


LA 


Antenna diameter 


3.7 m 


12.8 m 


Number of antennas 


10 




Baseline lengths (current) 


5 -20 m 


18 - 110m 


Primary beam at 15.7 GHz 


20'.1 


5' .5 


Synthesized beam 


a 3' 


rj 30" 


Flux sensitivity 


30 mJy s~ 1/2 


3 mJy s~ 1/2 


Observing frequency 


13.9- 18.2GHz 


13.9- 18.2GHz 


Bandwidth 


4.3 GHz 


4.3 GHz 


Number of channels 


6 


6 


Channel bandwidth 


0.72 GHz 


0.72 GHz 



Table 2. Assumed I+Q flux densities of 3C286 and 3C48, and errors on flux 
measurements in each frequency channel, over the commonly-used AMI 
SA bandwidth. 



Channel 


v/GHz 


5 3C2S6 /Jy 


s 3C48 /Jy 


<rs 


3 


14.2 


3.61 


1.73 


6.5% 


4 


15.0 


3.49 


1.65 


5.0% 


5 


15.7 


3.37 


1.57 


4.0% 


6 


16.4 


3.26 


1.49 


3.5% 


7 


17.1 


3.16 


1.43 


4.0% 


8 


17.9 


3.06 


1.37 


7.0% 



2 THE TELESCOPE 

AMI comprises two arrays: the Small Array (SA) which con- 
sists of ten 3.7-m diameter antennas, and the Large Array (LA) 
with eight 13-m antennas, located a t Lord's Bridge, Cambridge 
(AMI Consortiu m": Zwart et alj|2008l) . The higher resolution and 
flux sensitivity of the LA allows contaminating radio sources to 
be dealt with. These sources can then be subtracted from the 
SA maps. A summary of the technical details of AMI is given 
in Table ffl Further details on the telescope can be found in 
I A MI Consortium: Zwart et alj d2008h . 



3 OBSERVATIONS AND DATA REDUCTION 

Observations of Abell 2146 were made by the SA and LA be- 
tween November 2009 and March 2010, yielding approximately 
9 hours of good quality SA data; approximately the same amount 
of data suffered from artifacts and was discarded. Data reduction 
was perfomed using reduce, a local software too l developed for 
the Very Small Array (VSA) dWatson etal]|2003h and AMI (see 
e.g. lAMI Consortium: Zwart et a l. (2008) for further details). This 
package is designed to apply path delay corrections and a series of 
algorithms tailored to remove automatically bad data points arising 
from interference, shadowing, hardware and other errors. We ap- 
ply amplitude clips at a 3cr level. Periods where the data has been 
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contaminated by interference are excised. These interference sig- 
nals are identified as persistent high amplitude signals in the lag 
domain, which appear in all the lag channels. The system temper- 
ature is monitored by a modulated noise signal sent to the front- 
end of each antenna and synchronously detected at the end of each 
intermediate-frequency channel and is used in reduce to correct 
the amplitude scale on an antenna basis. If the system temperature 
falls below 10% of the nominal value of an antenna the associated 
datapoints a re removed. Fo r furth er details on the AMI reduction 
pipeline see iHurlev-Walken [2009). Additional manual flagging of 
remaining bad data points is done to ensure the quality of the data. 
The correlator data are then Fourier transformed into the frequency 
domain and stored on disk in FITS format. 

Flux calibration was performed using short observations of 
primary calibrators, either 3C48 or 3C286. The flux densities for 
3C48 and 3C286, see T abled are in agre ement with lBaars et al.l 
(1977) at 16-GHz. Since lBaars et alJdl977h measure I, as opposed 
to AMI which measures I + Q, the flux densities were corrected by 
interpolating from VLA 5-, 8- and 22- GHz observations. Previous 
tests have shown this calibration to be acc urate to better than 5 per 
cent ( I AMI Consortium: Scaife et al . 2009). The phase is calibrated 
using in terleaved calibrators selected from the Jodrell Bank VLA 
Surve y jPatnaik et al] 1 1992k iBrowne et all 1 19981: IWilkinson etai] 
1 19981) based on their proximity and flux density. The phase cali- 
brators used for the observations of Abell 2146 were J1642+6856 
for the SA and J 1623+6624 for the LA. These phase calibrators 
were interleaved approximately every hour for the SA and every 
ten minutes for the LA. 



3.1 Source subtraction 

Contamination from radio point sources at « 15 GHz can signifi- 
cantly obscure the SZ signal and must therefore be taken into ac- 
count in any SZ effect analyses at these frequencies. The higher 
resolution and flux sensitivity of the LA is exploited to deter- 
mine the position of the sources in the SA maps accurately in a 
short amount of time. Local maxima on the continuum LA maps 
above 4<x„, where cr„ is the corresponding value in Janskys per 
beam at that pixel in the noise map, are identified as LA de- 
tecte d sources using AMI-developed sourc e extraction software 
(see lAMI Consortium: Franzen et all ( l2010h ). Out of these LA- 
detected sources only those which appear within 0.1 of the SA 
power primary beam having an apparent flux above 4cr„ on the SA 
map are included in the source model. 

Every source in the source model is parameterized by a 
position, a spectral index and a flux density whose priors are 
based on the LA measurements. The source model is analysed 
by McAdam (Monte Carlo Astronomical Detection and Measure- 
ment), a Bayesian analysis pac kage for cluster detec tion and pa- 
rameter extr action developed by [Marshall et al. (2003) and adapted 
for AMI by iFeroz et al] d2009l) . which fits a probability distribu- 
tion to the source flux densities at the positions given by the LA. 
The source flux densities are fitted by McAdam to allow for pos- 
sible intercalibration difference between the two AMI arrays and 
for source variability. The mean source flux-density values are then 
used to subtract the sources from the SA map. 



4 BAYESIAN ANALYSIS OF CLUSTERS 
4.1 Bayesian inference 

The cluster analys is software implemented in this paper 
dMarshall et 111 12003) is based on Bayesian inference. This robust 
methodology constrains a set of parameters, 0, given a model or 
hypothesis, H and the corresponding data, D, using Bayes' theo- 
rem: 



Pr(0|D, H ) = 



Pr(D|0, H ) Pr(0|#) 
Pt(D\H) ' 



(1) 



Here Pr(0|D, H) = P(0) is the posterior probability distribution of 
the parameters, Pr(Z>|0,#) = £(Q) is the likelihood, Pr(0|#) = 
7r(0) is the prior probability distribution and Fr(D\H) = Z. the 
Bayesian evidence. If chosen wisely, incorporating the prior knowl- 
edge into the analysis reduces the amount of parameter space to be 
sampled and allows meaningful model selection. Bayesian infer- 
ence can serve as a tool for two main purposes: 

(i) Parameter estimation — In this case, the evidence factor can 
be neglected since it is independent of the model parameters, 0. 
Sampling techniques can then be used to explore the unnormalized 
posterior distributions. One obtains a set of samples from the pa- 
rameter space distributed according to the posterior. Constraints on 
individual parameters can then be obtained by marginalising over 
the other parameters. 

(ii) Model selection — The evidence is crucial for ranking mod- 
els for the data. It is defined as the factor required for normalising 
the posterior over 0: 



f £(®)n(0)d D 0, 



(2) 



where D is the dimensionality of the parameter space. This factor 
represents an average of the likelihood over the prior and will there- 
fore favour models with high likelihood values throughout the en- 
tirety of parameter space. This satisfies Occam's razor which states 
that the evidence will be larger for simple models with compact 
parameter spaces than for more complex ones, unless the latter fit 
the data significantly better. Deciding which of two models, Hq and 
Hi , best fits the data can be done by computing the ratio 



Pr(ffi|g) = PrfPIHQPr^,) = Za Piffl) 
Pi(H \D) Pr(Z>Iff )Pr(ff ) Z Pr(#o)' 



(3) 



where |^ is the prior probability ratio set before any conclusions 
have been drawn from the dataset. 



4.2 Nested sampling 

Nested sampling i s a Monte Carlo method introduced by 
ISkilling et al] d2004l) which focuses on the efficient calculation of 
evidences and generate s p osterior distributions as a by-pr oduct. 
IFeroz & Hobsonl d2008h and lFeroz. Hobson & Bridges! d2009l) have 
developed this sampling framework and implemented the Multi- 
Nest algorithm. This algorithm can sample from posterior distribu- 
tions where multiple modes and/or large (curving) degeneracies are 
present. This robust technique has reduced by a factor of as 100 the 
computational costs incurred during Bayesian parameter estima- 
tion and model selection. For this reason the analysis in this paper 
is based on this technique. 
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5 PHYSICAL MODEL AND ASSUMPTIONS 
5.1 Interferometric data model 

An interferometer, like AMI, operating at a frequency, v, measure 
samples from the complex visibility plane / v («). These are given by 
a weighted Fourier transform of the surface brightness, I v (x): 



/y(«) = J 



A v (x)Iy(x)exp(2niu ■ x)d x, 



(4) 



where x is the position relative to the phase centre, A v {x) is the 
(power) primary beam of the antennas at an observing frequency, v 
(normalized to unity at its peak) and u is the baseline vector in units 
of wavelength. In our model we assume the measured visibilities 
can be defined as 



V V (U,) = /,(«,) + Ny(Ui), 



(5) 



where I v (u) is the signal component, which contains contributions 
from the cluster SZ effect signal and identified radio point sources 
and 2Vy(»,-) is a generalized noise component that includes signals 
from unresolved point sources, primordial CMB anisotropics and 
instrumental noise. 



5.2 Cluster models 

In order to calculate the contribution of the cluster SZ signal to the 
visibility data the Comptonization par ameter of the clus ter, y(s), 
across the sky must be determined (see lFeroz et al. I j2009h for fur- 
ther details). This parameter is the integral of the gas pressure along 
the line of sight / through the cluster: 



y(s) 



—^j I n e k B Tdl<x f pgTdl, 



(6) 



where crj is the Thomson scattering cross-section, n c is the electron 
number density, which is derived from equation 10, m c is the elec- 
tron mass, c is the speed of light and k% is the Boltzmann constant. 
s = 6D g is the deprojected radius such that r 2 = s 2 + I 2 and D s is 
the angular diameter distance to the cluster which can be calculated 
for clusters at redshifts, z, using 



(7) 



(1+z) 

We set r ]im in equation ((6} to 20/i~'Mpc — this result has been 
tested and shown to be large enough even for small values of /3 
dMarshall et aT]|2003l) . 

The cluster geometry, as well as two linearly independent 
functions of its temperature and density profiles, must be specified 
to compute the Comptonization parameter. For the cluster geome- 
try we have chosen a spherical cluster model as a first approxima- 
tion. The temperature profile is assumed to be constant throughout 
the cluster. An isothermal /?-model is assumed for the cluster gas 
density, p g dCavaliere & Fusco-Fermianoll 19781) : 

Pg(0) 



i + (t) 



21 T 



where 



p g (r) = /U e tte(r), 



(8) 



(9) 



p c = 1 . 14ra p is the gas mass per electron and m p is the proton mass. 
The core radius, r c , gives the density profile a flat top at low j- and 
p g has a logarithmic slope of 3/3 at large — . 



Parameter estimates can depend on the way the cluster model 
is parameterized. We examine the impact of different physical as- 
sumptions by presenting the parameter estimates for Abell 2146 
obtained using three different cluster parameterizations (or 'mod- 
els'). Modelled sources for all three models are characterised by 
three parameters: position, flux density and spectral index. The cor- 
responding priors for these parameters are given in Section [5. 3. 21 
The parameterizations of the sources and the source priors are the 
same in all three models, unlike the cluster parameterizations which 
do change between models. The mean values fitted by our McAdam 
software to both the source and cluster sampling parameters will, 
however, vary for each cluster model. We proceed to describe our 
three cluster parameterizations and their results. 

Tables [4] and [3] indicate which parameters are derived in each 
model and the assumptions made in each case. A summary of the 
sampling parameters for each model together with their priors is 
given in Table [5] 

5.2.1 Cluster model 1 

Our first model, henceforth Ml, is based on traditional methods for 
the analysis of SZ and X-ray data. The sampling parameters for Ml 
are: 

• (*c, ^c) — the position of the cluster centroid on the sky. 

• T — the temperature of the cluster gas, which is assumed to be 
uniform. 

• p — defines the outer logarithmic slope of the /^-profile. 

• r c — gives the density profile a flat top at low r 

• M g (r 2 oo) — the gas mass inside a radius, r 2 oo, which is the ra- 
dius at which the average total density is 200 times p cr j t , the critical 
density for closure of the Universe. 

• z — the cluster redshift. 

In applying this cluster model to Abell 2146 both z and T are 
assumed to be known, which is equivalent to assigning them delta- 
function priors (see Table[5j ■ 

The derived parameters for Ml are: 

• r x — the radius at which the average total density is X times 

Pcrif 

• Mx(rx) — the total cluster mass within the radius r x . 

• M g (r x ) — the cluster gas mass within the radius r x . 

• fg( r x) — the cluster average gas fraction within the radius r x . 

• p g (0) — the central gas density. 

• y(0) — the central Comptonization parameter. 

In this model, the cluster gas is assumed to be in hydrostatic 
equilibrium with the total gravitational potential of the cluster, O, 
which is dominated by dark matter. As a result, the gravitational 
potential must satisfy 

d<5 _ 1 dp 

dr p g dr 



(10) 



This equation can be simplified if the cluster gas consists purely of 
ideal gas with a uniform temperature, T, to give 

dlogp„ Gp M T (r) 



dlogr 



k B T 



(11) 



where n is the ma ss per particle, p 



0.6ra n 



rnPe (see 



iMarshall et afll2 003). Expressions for the total mass of the cluster, 
Mj(r x ), can be obtained for spherical symmetry 

4k 



M T (r x ) : 



~r X XPcrili 



(12) 
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or by integrating the isothermal /3-model for the density profile in 
CCD, 



MtOx) 



4 3pk B T 



(13) 



Combining equations l l 1 2t and d 1 3 1 > leads to an expression for r x , 



4^GXp cri , 



(14) 



The total mass of the cluster within a certain radius, Mj(r x ), is sub- 
sequently determined by substituting r x into equation d!2t . Once 
Mj(r x ) and M g (r x ) are known, the gas fraction, f g (r x ), can be com- 
puted using the relation 



/«(*) : 



M g (r x ) 



(15) 



M T (r x ) 

We consider values for X = 200 and 500. For X = 500, M g (r x ) is 
not a sampling parameter but is calculated using the expression 

C rx Anr' 1 

M g (r x )=p g (0) w dr', (16) 

Jo r ,,"|t 

[ 1+ f] 

Also, Pu(0), in equation (8), can be recovered by numerically inte- 
grating the gas density profile up to r 20 o, equation d 1 6b . and setting 
the result equal to M g (r 2 oo). 



J. 2. 2 Cluster model 2 

Our second model, M2, has the same sampling parameters as 
Ml with the exception of T, which becomes a derived parameter, 
and / g (r 2 oo), which becomes a sampling parameter. Sampling from 
/ g (r 2 oo) and M g (r 20 o) allows M T (r 20 o) to be calculated using equa- 
tion d 1 5 b - r 2 oo can then be computed simply by rearranging equation 
d 1 2b . The temperature of the cluster gas can be obtained by com- 
bining equations i ll It and <[Sj to yield 



GyU M T (r 200 )(r 2 c + r 2 2Q0 ) 
3k B /3 r\ m 



(17) 



which is based upon the assumption that the cluster is in hydrostatic 
equilibrium and described well by a /3-profile. The derived parame- 
ters at r 500 are calculated in the same way as in Ml ; once M g (r 500 ) is 
obtained from equation \16\ and rsoo from equation d!4t , Mt^oo) 
is calculated by assuming the cluster is spherical, equation U2b . 
fg( r 50o) can then be recovered using the relation in dl5l > ■ 

5.2.3 Cluster model 3 

In the third model, M3, the sampling and derived parameters are 
the same as in M2. The only difference between M2 and M3 is the 
way T is calculated. M3 uses an M-T relation to derive T which 
allows T to be obtained without relying on the cluster being in hy- 
drostatic equilibrium, a necessary assumption in M2. Moreover, at 
r 2 oo, all the other cluster parameter estimates of M3 are free from 
the assumption of hydrostatic equilibrium. However, this assump- 
tion needs to be made to obtain cluster parameters at 7-500 (see Sec- 
tion |5.2.2t . 

If the cluster is assumed to be virialized and to contain a small 
amount of unseen energy density in the form of turbulence, bulk 
motions or magnetic fields, the average cluster gas temperature, T, 
can be obtained using the mass-temperature (M-T) relation for a 
singular, isothermal sphere (SIS) based on the virial theorem, 



G/jM t 

2r 20 o 



! (^(200p cri ,)) 



1/3 M T 



2/3 



: 8.2keV 



Mi 



10 I5 /i-'M o / \ H ) 



2/3 



lH(z)\ 2 ' 3 



(18) 
(19) 

(20) 



where H is the Hubble parameter. In our cluster model we use the 
well-behaved /J-profile, equation l[8j, rather than the SIS density 
profile which is singular at r = 0. This different choice for the den- 
sity profile will introduce a factor to the M-T relation in equation 
d20b . From cluster simulations we find that this factor varies be- 
tween 0.7 - 1.2. 



5.2.4 M-T relation and hydrostatic equilibrium 

The results obtained from running McAdam with three different 
models are useful for assessing the validity of some of the assump- 
tions made in each model. Traditional models tend to assume clus- 
ters are isothermal, spherical, virialized and in hydrostatic equi- 
librium. All of these assumptions are particularly inappropriate for 
cluster mergers like Abell 2146. The first two assumptions are made 
in the three models presented in this paper to simplify the cluster 
model; but note that the spherical assumption is not bad here be- 
cause our SZ measurements are sensitive to the larger scales of the 
cluster. 

M2 also assumes hydrostatic equilibrium to obtain an esti- 
mate for T. After the gravitational collapse of a cluster, the hot 
gas in the ICM tends to reach equilibrium when the force exerted 
by the thermal pressure gradient of the ICM balances that from the 
cluster's own gravitational force. An underlying assumption is that 
the gas pressure is provided entirely by thermal pressure. In real- 
ity, there are many non-thermal sources of pressure support present 
in most clusters such as turbulent gas motions which can provide 
~ 10 - 20% of the total pressure support even in relaxed clusters 
dSchuecker. Bohringer and Vogedl2004l: iRasia et alj|2006h . In the 
case of Abell 2146, a complex merging syste m with two detecte d 
shocks propagating at ss 1900 and 2200kms _1 dRussell et al.l2010h . 
there is significant non-thermal pressure support provided by bulk 
motions in the ICM. 

Relating radius, temperature and total mass via the virial the- 
orem in practice also assumes that the kinetic energy is in the form 
of internal energy of the particles, as evidenced by the SZ signal, 
so that turbulent motions, bulk motions and everything else are ig- 
nored. But this use of the virial theorem has an advantage over hy- 
drostatic equilibirum in the case of Abell 2146 since our knowledge 
of the mass ratio of the two merging systems enables us to set a 

limit on the degree to which the use of the M-T relation, T cc M^ , 
bias es our temperature es timate. 

iRussell et alJ j2010h find the fractional mass of the merging 
cluster to be between 25 and 33 percent, in which case the average 
temperature of the merging system will be as 10% higher when all 
the gas mass of the subcluster has merged with that of the primary 
cluster than prior to the start of the merger event. Therefore, pro- 
vided the primary cluster was virialized pre-merger, our estimate 
for T using the M-T relation in equation {20} is little affected by 
the merger. 
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Table 4. Summary of the main assumptions made in the calculation of the derived parameters for each model. H stands for hydrostatic equilibrium, M-T for 
the mass-temperature relation given in equation J20t , B for isothermal /3-profile, S for spherical geometry and N/A means not applicable, since that parameter 
is a sampling parameter for that particular model. 



Model Assumptions 


Derived Parameter 


Model 1 


Model 2 


Model 3 


r 2 Qo/h-' Mpc 


H, S, B; Eq.rm 


S;Eq[l2] 


S; EqfHl 


M T (r 20 o) 


S;Eq[l2] 


Eq[B] 


Eq[l5] 




Eq[B] 


N/A 


N/A 


M g (r 50 o)//r 2 M 


S, B; Eq[T6l 


S.B:Eq[T6l 


S, B; Eq[T6l 


r 5Q0 /h- 1 Mpc 


H, S, B; Eq[T4l 


H, S, B; EqfTil 


H, S, B; Eq[T4l 


M T (r 50 o) 


S;Eq[l2] 


S;Eq[l2] 


S;Eq[l2] 


f g (r 50Q )/h- [ 


Eq[l5] 


Eq[B] 


Eq[T5] 


TkeV 


N/A 


H;Eq[T7l 


M-T; Eq[20l 



Table 3. Summary of the derived parameters for each cluster model. 



Derived Parameter 


Model 


loo & rsoo/h- 1 Mpc 


All 


M T (loo) & M T (r5oo)//i _1 M G 


All 


M g 0- 500 )//r 2 M o 


All 


y All 


«e 


All 


TkeV 


M2, M3 


/gdooV/r 1 


Ml 


./gteoo)//^ 1 


All 



5.3.2 Source priors 

As with the cluster priors, the source priors are assumed to be sep- 
arable, such that 

7! (0 S ) = n (x s ) n (v s ) 7i (So) n (a) . 

7r(x s ) and /r(y s ) are given delta priors at the source position 
found from the high-resolution LA maps. The flux-density priors 
for modelled sources on the other hand are chosen to be Gaussians 
centered on the flux-density value given by the LA with a as 40% 
of the LA source flux. Tight constraints on the flux-density pri- 
ors are best avoided due to inter-array calibration differences and 
source variability. The channel flux densities taken from the LA 
data are used to calculate an estimate for the spectral index of each 
source. The spectral index prior is then set as a Gaussian centered 
at the predicted LA value with a width cr = 1 . 



5.3 Priors 

5.3.1 Cluster priors 

For simplicity the priors are assumed to be separable. The priors 
used in the analysis of Abell 2146 are given in Table H] 

We note that, although the prior on M g (r 2 oo) assumes the clus- 
ter produces a non-zero SZ effect, it is wide enough that our results 
will not be biased. In fact, our posterior distributions for M g (r 50 o) 
peak at M g (r 50 o) > 3 x 10 I3 //T 2 M G and have fallen to zero by 
M g (r 500 ) > 2 x 10 l3 //i~ 2 M o , while our prior for M g (r 2 oo) extends 
down to 1 x 10 13 //^ 2 M o 

The prior on the gas mass fraction was set to a Gaussian 
centered at the WMAP7 best-fit value, / g = 0.12ft -1 , with cr = 
0.016/T 1 . This result was obtained from WMAP7 estimates of 
= 0.266 + 0.029, Cl b = 0.0449 ± 0.0028 and h = 0.710 ± 0.025 
using the relation f b = where f b is the universal baryon frac- 
tion (|LMSonet^][20l3). The prior on f g can be based on f b since 
/ g in clusters at large radii approaches f b . The prior on the position 
of the cluster was a Gaussian with a = 60" centered at the X-ray 
centroid. 



6 RESULTS 

6.1 Maps and evidences 

Fifteen sources were detected above 4cr„ on the LA map, Fig. [T] 
McAdam was used to determine the flux densities and spectral in- 
dices of these sources in the S A data. The standard AIPS tasks were 
used to clean the images with a single clean box. No primary beam 
correction has been applied to the AMI maps presented in this paper 
such that the thermal noise, o~ n , is constant throughout the map. The 
task imean was applied to the data to determine the noise level on 
the maps. Contours increasing linearly in units of cr„ were used to 
produce all the contour maps. The half-power contour of the syn- 
thesized beam for each map is shown at the bottom left of each 
map. 

Further analysis was undertaken in the visibility plane taking 
into account receiver noise, radio sources, and contributions from 
primary CMB imprints. Figs [2] and [3] show the SA maps of Abell 
2146 before and after radio source subtraction. The source sub- 
traction was performed at the LA source position using the mean 
flux-density estimates given by the McAdam results of M3, Table[6] 
Sources with a high signal-to-noise ratio and close to the pointing 
centre tend to have good agreement between flux densities mea- 
sured by the LA and those obtained by McAdam. Possible reasons 



Sunyaev-Zel'dovich observation ofAbell 2146 7 

Table 5. Summary of the priors for the sampling parameters in each model. 



Parameter 


Models 


Prior Type 


Values 


Origin 


*c,;yc" 


all 


gaussian at *x-ray> 0" = 60 


15 h 56 m 07 s , +66°2i:35" 


Ebeline et al. (2000) 


P 


all 


uniform 


0.3-2.5 


Marshall et al. (2003) 


M g (r 200 )//r 2 M G 


all 


uniform in log 


10 13 - 10 15 


physically reasonable 


r c //i~'kpc 


all 


uniform 


10 - 1000 


physically reasonable 


z 


all 


delta 


0.23 


Ebelins et al. (2000) 


/g(r2oo)/T 1 


M2, M3 


gaussian, cr = 0.016 


0.12 


Larson et al. (2010) 


TkeV 


Ml 


delta 


6.7 


Russell et al. (2010) 



for source flux-density discrepancies between the arrays, in partic- 
ular for the remaining sources, include: a poorer fit of the Gaussian 
modelled primary beam at large m'-distances from the pointing cen- 
tre, loss of signal due to the white light fringes falling off the end 
of the correlator, time and bandwidth smearing, correlator artifacts, 
source variability and, some sources with low signal-to-noise ratios 
detected on the LA, might appear as noise features on the SA. 

It should be noted that, since a single flux-density value is 
used for subtracting the modelled sources in the map-plane, the 
radio source subtracted map does not reflect the uncertainty in 
the McAdam derived flux-density estimates Q Nevert heless, flux- 
densit y estimates given by McAdam have been tested in lFeroz et al.l 
J2009h and shown to be reliable. Fig. [9] shows that there is no de- 
generacy between the flux density fitted for source A in Fig.[3]and 
the fitted values for M g (r2oo), the cluster gas mass within r2oo The 
detection of Abell 2146 in the AMI data is confirmed by compar- 
ing the evidence obtained by running McAdam with a model in- 
cluding SZ + CMB primordial structure + radio sources + receiver 
noise and the null evidence, which corresponds to a model with- 
out a cluster, i.e. simply CMB + radio sources + receiver noise. 
The first model, which included an SZ feature, was found to be e 15 
times more probable than one without. 

In Fig. [4] a 0.6-kyi taper is used to enhance large scale struc- 
ture and consequently the signal-to-noise ratio of the SZ effect. The 
peak decrement in it is as 13cr. 

The AMI SZ maps are compared to the Chandra X-ray emis- 
sion and projected temperature maps for Abell 2146 in the discus- 
sion, Section|7j 



6.2 Parameter estimates from three cluster models 

McAdam was run on the same Abell 2146 data for each of the three 
models described in |5.2| The results obtained for these models are 
shown in Figs [5] to [TTJ The contours in all the 2D marginalized 
posterior distributions represent 68% and 95% confidence limits. 
Axis labels for M g (r 2 oo) are in units of 10 13 for clarity. 



Note that, unlike for the radio source subtracted maps, when obtaining 
estimates for the cluster parameters the whole probability distribution for 
the source flux density is taken into account, such that a larger uncertainty 
in the source flux densities will lead to wider distributions in the cluster 
parameters. 



66 45 




1 


i i i i 




40 

35 






TO, 




30 






' 3 
>« 




DECLINATION (J2000; 






a . S 




10 










05 
00 










65 5! 


m 1 


i i i 


i i i i 


i- 



16 00 15 59 58 51 56 55 54 53 52 

EIGHT ASCENSION (J2000J 
Cont peak flux = 5.8479E-03 JY/BEAM 
Leus = 1.806E-04 * f-10, -9, -8, -7, -6, -5, -4, 
-3, 3, 4, 5, S, 7, 8, 9, 10) 



Figure 1. LA contour map. 




Figure 2. SA map before source subtraction. The crosses indicate the posi- 
tion of the sources detected on the LA map. 
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Table 6. List of the detected sources with their J2000 position coordinates, as determined by the LA map. Columns 3 and 4 show the flux densities of the 
detected sources at 16-GHz (S a) given by McAdam using M3 with their associated Gaussian errors. For comparison, the LA measured flux densities at the 
same frequency are given. The letters represent the labeled sources in Fig. [3] 



Source 


RA (h m s) 


Dec (" ' ") 


McADAM-fitted S i6(mjy) 




LA S 16 (mJy) 


A 


15 56 04.23 


+66 22 12.94 


5.92 


0.18 


5.95 


B 


15 54 30.95 


+66 36 39.58 


0.60 


0.29 


0.61 


C 


15 56 14.30 


+66 20 53.45 


1.83 


0.14 


1.70 


D 


15 56 36.51 


+66 35 21.65 


2.15 


0.15 


1.65 


E 


15 55 57.42 


+66 20 03.11 


1.65 


0.08 


1.64 


F 


15 58 10.23 


+66 24 35.72 


1.49 


0.12 


1.29 


G 


15 54 03.96 


+66 28 41.90 


1.12 


0.15 


0.74 


H 


15 55 25.67 


+66 22 03.96 


0.48 


0.05 


0.67 


I 


15 55 10.84 


+66 19 45.82 


0.61 


0.06 


0.65 


J 


15 57 09.46 


+66 22 37.62 


0.43 


0.06 


0.63 


K 


15 54 47.50 


+66 28 37.43 


0.91 


0.09 


0.53 


L 


15 54 49.11 


+66 14 21.49 


0.72 


0.09 


0.47 


M 


15 56 15.40 


+66 22 44.48 


0.16 


0.07 


0.43 


N 


15 56 27.90 


+66 19 43.82 


0.11 


0.05 


0.33 


O 


15 57 56.10 


+66 22 49.80 


0.30 


0.07 


0.49 
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RIGHT ASCENSION (J2000) 
Cont peak flux = -1.3725E-03 JY/BEAM 
Levs = 1.40BE-04 * (-10, -B, -6, -4, -2, 2, 4, 6, 

B , 10, 12) 



Figure 3. SA contour map after source subtraction. The letters represent the 
position of the sources detected on the LA map. 
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Cont peak flux = -1.9027E-03 JY/BEAM 
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-3, -2, 2, 3, 4, 5, 6, 7, B, 9, 10) 



Figure 4. SA map after source subtraction using a 0.6-t-l taper. The crosses 
represent the position of the sources detected on the LA map. 



6.2.1 Cluster model 1 

The 2-D and 1-D marginalized posterior probability distributions 
for the parameters of Ml are depicted in Figs.[5]and[6] respectively. 

Ml is representative of the more conventional method for ex- 
tracting cluster parameters from SZ data. In this model, the aver- 
age cluster gas temperature within r 2 oo is assumed to be known, 
from X-ray measurements, allowing the morphology of the cluster, 
namely r x , to be inferred by assuming the cluster is spherical, in hy- 



drostatic equilibrium and described well by an isothermal y3-model. 
The overall bias on r x arises from all of these assumptions, which 
are particularly unphysical in a cluster merger like Abell 2146, and 
is therefore expected to be large. Indeed, by comparing Figs.|6]and 
II llwe find that r 2 m is overestimated with respect to the value ob- 
tained in M3, our most physically motivated model. Moreover, in 
Ml, M g (r x ) for X = 500 and 1000 depends on r x , which results in 
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Figure 5. Two-dimensional marginalized posterior distributions for the 
sampling parameters of Abell 2146 — Ml. 
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Figure 6. One-dimensional posterior probability distributions for selected 
derived parameters of Abell 2146 — Ml. We note that the axes for the plots 
of the 1-D and 2-D marginalized posterior distributions of both the sampling 
and derived parameters are tailored to suit the results of each model and will 
therefore be different in each case 



the bias on r x to be propagated to the remaining derived parameters 
for these values of X. 



6.2.2 Cluster model 2 

The 2-D and 1-D marginalized posterior probability distributions 
for the parameters of M2 are depicted in Figs|7]and[8] 

M2 introduces a new sampling parameter, / e (r 2 oo)- Sampling 
from this parameter allows more prior information to be included 




Figure 7. Two-dimensional marginalized posterior distributions for the 
sampling parameters of Abell 2146 — M2. 
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Figure 8. One-dimensional marginalized posterior distributions for the de- 
rived parameters of Abell 2146 — M2. 



in the analysis, which has the effect of constraining the parame- 
ter distributions better than in Ml. It has a great advantage over 
Ml, namely, the only parameter obtained by assuming hydrostatic 
equilibrium is the temperature, which is not used explicitly in the 
calculation of the other derived parameters at r 2 oo- 

Fig. [9] shows the two-dimensional marginalized posterior dis- 
tribution for the flux density of source A, Sa, and M g (>2oo)-we 
choose to plot 5a since source A is the brightest source close to the 
pointing centre. One can see from Fig. [9] that S a and M a (r 2 oo) do 
not appear to be significantly correlated. This is confirmed by the 
sample correlation, which was found to be 0.12. We note that the 
sample correlation remains unaffected by shifts of origin or changes 
of scale in Sa and M g (>2oo)- The flux density of source A is given 
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Figure 9. Two-dimensional marginalized posterior distribution for the flux 
of source A shown in Fig. [3] Sa, and the cluster gas mass within r2oo, 
M g (r 2 oo). 

a Gaussian prior and yet the LA-measured and McAdam-derived 
flux-density estimates for this source are very close. 

6.2.3 Cluster model 3 

The 1-D and 2-D marginalized posterior probability distributions 
for the parameters of M3 are presented in FigsllOlandfTTI 

The only difference between M2 and M3 is in how the average 
cluster gas temperature at r 2 oo, T, is calculated. To obtain an esti- 
mate for T, M2 assumes the cluster is in hydrostatic equilibrium 
while M3 uses the M-T relation in equation l !20t . which assumes 
the cluster is virialized and contains no unseen energy density. 



7 DISCUSSION 

7.1 Comparison with X-ray maps 

Two new Chandra observ ations of Abell 2146 were taken in April 
2009 dRussell et alj20l6T) . Fig.|13|shows the exposure-corrected X- 
ray image taken in the 0.3 - 5.0 keV energy band smoothed with 
a 2D Gaussian of cr = 1.5 arcseconds superimposed with the AMI 
SZ effect from Fig. [3] The AMI av-coverage is well-filled and goes 
down to « 180/1 which corresponds to a maximum angular scale of 
ss 10 arcminutes or a cluster radius of as 1.1 Mpc. Thus, in practice, 
the SZ signal traces a more extended region of the gas than the 
X-ray data. Any small features in the cluster environment are not 
resolved by the SA maps which consequently appear much more 
uniform than the X-ray maps. Nevertheless, given the synthesized 
beams in Figs.|3]and|4]the SZ effects in these two figures appear to 
show signs of some real extended emission. To verify that we have 
resolved the SZ decrement we bin the data from the CLEANed, radio 
source subtracted, non-tapered map of Abell 2146, Fig. [2] in bins 
of 100/i and plot it against baseline, see Fig.[T2] The signal steadily 
becomes more negative from scales of 800/1 to 200/1; it is on these 
larger scales that we find the most negative binned value for the SZ 
decrement, demonstrating the sensitivity of the SA to large angular 




Figure 10. Two-dimensional marginalized posterior distributions for the 
sampling parameters of Abell 2146 — M3. 
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Figure 11. One-dimensional marginalized posterior distributions for the de- 
rived parameters of Abell 2146 — M3. 



scales. To determine the shape of the cluster in greater detail high 
resolution SZ observations are needed. 

During a cluster merger, elongations in the dark matter and 
gas components are expected. In general, the orientation of this 
elongation for both components tends to be parallel to the merger 
axis, though the gas component can also be extended in a direction 
perpendicular to the mer ger ax is due to adiabatic compressions in 
the ICM (Roettiger et al. 1997), as shown in simulations of cluster 
mergers dPoole et alj|2007h . We fitted a six-component (position, 
peak intensity, major and minor axes and position angle) elliptical 
Gaussian to the SZ decrement in our 0.6 tapered map, Fig. [4] and 
a zero level using the AIPS task JMFIT. The results for the param- 
eters defining the shape of the fitted ellipse are given in Table [7] 
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Figure 12. Binned data from the CLEANed, radio source subtracted, non- 
tapered map of Abell 2146 in bins of 100/1 against baseline (in kA). It 
should be noted that the FWHM of the aperture illumination function of 
the AMI SAis a 185 A such that the visibilities in each bin are not entirely 
independent. The baseline distance corresponding to the McAdam derived 
parameters rioo, '"soo and nooo were found to be 0.421oi,0.816k/l, 1.46k/l, 
respectively. 



Table 7. IMFIT results for the parameters of the ellipse fitted to the SZ 
decrement in the 0.6-tapered SA CLEANed maps. The extension of the minor 
and major axes are given in arcseconds and the position angle in degrees. 





Nominal 


Minimum 


Maximum 


Major axis 


205 


171 


236 


Minor axis 


145 


109 


175 


Position angle 


46 


3 


68 



The nominal results indicate that the semi-major axis has a posi- 
tion angle of 46°. The orientation of the SZ signal along this axis 
seems to be ^orthogonal to the elongation of the X-ray signal, see 
Fig. [13] Shock fronts like the ones observed in Abell 2146 can only 
be detected during the early stages of the merger, before they have 
reached the outer regions of the system which suggests that the gas 
disturbances in the cluster periphery are less intense than those near 
the dense core. 

This is supported by the different signal distributions of the X- 
ray and SZ effect data. The gas is relatively undisturbed in the clus- 
ter periphery while in the inner regions the core pas sage has dis- 
place d the local gas at right angles to the merger axis dRussell et al.l 

[2oToh . 

The total mass can also be estim ated from the X-ray M T (r 50 o) _ 
T relation (eg. VikMi nin et al.l d2006j) ) (note that here we use a dif- 
ferent scaling relation than elsewhere since we are concerned with 
cluster parameters at r 50 o). Excluding the cool core region, the X- 
ray spectroscopic temperature is 7.5 ± 0.3 keV, which corresponds 
to a mass M T (r 50 o) ~ 7±2x 10 I4 M o (using h 10 = 1.0). This method 
will likely overestimate the cluster mass as we expect the tempera- 
ture to have been temporarily boosted duri n g this major merger b y 
a factor of a few dRicker & SarazirJ fcOOlh . Rand all et alj d2002h ). 
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Figure 13. Chandra X-ray image superimposed with AMI SA SZ effect 
(no taper). The SA map is shown in black contours which go from —1.4 
mlybeam -1 to 0.001 mlybeam -1 in steps of +0.2 mJybeaiTT 1 . The grey 
scale shows the exposure-conected image in the 0.3 - 5.0 keV energy band 
smoothed by a 2D Gaussian cr = 1.5 arcsec (North is up and East is to the 
left). The logarithmic scale bar has units of photons cm~ 2 s~'arcsec~ 2 . 

A mass estimate for the Bullet cluster from the Mj(rsoo) - T re- 
lation produced a result approximately a fact or of 2.4 higher tha n 
the weak lensing result for the same region dMarkevitchl J2006)). 
If we assume the X-ray mass estimate for Abell 2146 is overes- 
timated by a similar factor, the cluster mass should be closer to 
M T (r 5 oo) ~ 3 x 10 14 /j~'M g , which is comparable with our SZ ef- 
fect result. However, simulations show that the transient increase in 
the X-ray temperature is dependent on the time since the collision, 
the impact para meter of the merger and th e mass ratio of the merg- 
ing clusters (e.g. lRitchie & T homas (2002)), which will be different 
for the Bullet cluster. A weak lensing analysis using new Subaru 
Suprime-Cam observations will produce a more accurate measure 
of the mass for comparison with the SZ effect result. 

7.2 Comparison with the 4.9-GHz VLA maps 

The VLA radio image taken at 4.9-GHz (NRAO/VLA Archive 
Survey) and the contours representing the LA map are superim- 
posed on the X-ray image in Fig. Q3] The presence of a bright 
source on top of the dense cluster core obscures any possible high- 
resolution SZ features in the LA map. High-resolution SZ images 
using the LA would be possible if higher resolution data taken at 
16-GHz were available for source subtraction. The longer baselines 
of the LA proved insufficient to remove the contaminant sources 
and no SZ effect decrement was seen on the source subtracted 
LA maps. High-resolution SZ effect measurements are necessary 
to disentangle the density and temperature distributions properly. 
These observations in o ther cluster mergers like the Bullet cluster 
dSawasachi et al. 2010) have revealed structure in the gas pressure 
distribution and are powerful tools for understanding the evolution 
of galaxy clusters. 

Radio halos are faint, large-scale sources that often span the 
entire cluster and are typically found in cluster mergers. Two hours 
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Right Ascension (J2000) 

Figure 14. Projected temperature map (keV) lRusselletai]|2010l) overlaid 
on black contours representing the SA SZ effect decrement. The contours 
go from -1.5 mJybeairT 1 to 0.1 mjybeairr' in steps of +0.1 mJybeairr' 
and the grey linear scale indicates the temperature variation in keV. 

of VLA observations in two configurations, C and D, towards 
A520 revealed a ra dio halo with a power of 6.4 x 10 24 W Hz -1 
dGovoni et alfeOOll) at 1.4-GHz. The Bullet cluster was also found 
to hav e a radio halo with a power of (4.3±0.3)x 10 25 WHz 1 at 1.3- 
GHz dLiang et al.l2000h . No low frequency radio data are currently 
available for Abell 2146. 4.9-GHz VLA observations of Abell 2146 
do not show signs for a radio halo, Fig. [15] though deeper observa- 
tions, particularly at lower frequencies where radio halo emission 
tends to be stronger, would be needed to determine whether a radio 
halo is present in Abell 2146. Sinc e such h a los ar e characterized 
by a s teeply falling spectrum (e.g iHanischl jl980h : iGovoni et al.l 
(2004)) and no radio halo emission was detected at 4.9-GHz, we 
do not expect our observations to be contaminated by this diffuse 
emission. 

A520 and 1E0657-56 are the only two clusters that have been 
found to have both bow shocks and radio halos. They have pro- 
vided unique information that allows determination of what pro- 
portion of the ultrarelativistic electrons producing the radio halo 
are generated as a res ult of merger-drive n turbulence, as o ppose d 
to shock acceleration jlvlarkevitch et alj|2002l ; lMarkevitchl l2006). 
Since Abell 2146 is the third cluster merger known to contain sub- 
stantially supersonic shock fronts, finding a radio halo would sig- 
nificantly improve our current understanding of how they are gen- 
erated and powered. 



7.3 Cluster Parameters 

The cluster parameters obtained from M3, our preferred model, are 
discussed below. 



7.3.1 Position 

The mean value for the position, RA 15 h 56 m 07 s Dec + 66°2i:33'.', 
with errors of 6 and 7 arcseconds respectively, coincides with the 
X-ray centroid position, RA 15 h 56 m 07 s Dec + 66°2i:35'.', as shown 
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Figure 15. VLA 4.9-GHz map in thick, grey contours overlaid on the 
AMI LA map, in thin, black contours and the X-ray grey map from Chan- 
dra observations. The logarithmic grey scale corresponds to the exposure- 
corrected X-ray image taken in the 0.3 - 5.0 keV energy band smoothed 
with a 2D Gaussian of a = 1.5 arcseconds and it is in units of photons 
cm~ 2 s~'arcsec~ 2 . The VLA and LA contours range from 0.5 mJybeairT 1 
to 9 mJy beam -1 in steps of 0.3 mJy beam -1 . 

in Fig.Q/3] However, the peak of the X-ray flux is significantly dis- 
placed from the peak of the SZ signal, as depicted in Fig.Q/3] The 
X-ray spectral luminosity is proportional to Jn 2 r~ 1/2 d/, while the 
SZ effect is a measure of the integrated line-of-sight pressure and 
is proportional to j n e Tdl. Therefore, the X-ray emission is more 
sensitive to substructure than the SZ data and peaks at the position 
of the dense cluster core. 



7.3.2 j3 and r c 

Results from running McAdam on large samples of clusters have 
forced the prior on p to be relaxed to include higher values (see 
lAMI Consortium: Zwart et alj d2008h ). The distributions for /? tend 
to favour higher values than typical X-ray estimates. However, this 
discrepancy is not surprising since previous studies have revealed 
incompatibilities in the f3 fits between X-ray and SZ effect profiles 
due to their diffe rent dependencies on parameters such as tempera- 
ture and density lHall man et al.l2007l) . The results show the degen- 
eracy between r c and /3 but also show evidence of strong constraints 
on this relation. This relation is positively correlated in Ml and M2, 
where the assumption of hydrostatic equilibrium is made to esti- 
mate parameters at r 2 oo, and negatively correlated in M3 where this 
assumption is avoided. 

7.3.3 Gas fraction 

Sampling from / g (f"2oo) allows further prior information to be intro- 
duced into the model which leads to better constrained parameter 
estimates. 

All the models were run through McAdam without data to 
check the effect of the priors on the results. From this test we dis- 
covered that in Ml the seemingly inconspicuous priors on the sam- 
pling parameters lead to an effective prior on / g that peaks around 
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0.01 and strongly disfavours values of f g « 0.1. Since in our current 
models SZ data alone cannot place strong constraints on /„, the ef- 
fective prior biases low the estimates of / s obtained in Ml. On the 
other hand, when running M3 without data, the effective prior on 
/ g does not change significantly from the Gaussian prior it was ini- 
tially given. Given the importance of analysing cluster models with- 
out any data to interpret their results, a detailed discussion of these 
no-data runs and the effects of cluster parameterization are pre- 
sented in the forthcoming paper, AM I Consortium : Qlamai e"et al.l 
42011b . 



7.3.4 Temperature 

The average cluster gas temperature within ^200 ~ 900 kpc for hjo = 
1.0 was found to be 4.5 + 0.5 keV. The projected emission- weighted 
temperature map, Fig. [14] shows a range of X-ray temperature mea- 
surements in different regions of the cluster. At the position of the 
most negative value of the SZ decrement, the X-ray temperature is 
ss 8 keV where as at a radius of ^ 500 kpc the temperature drops 
below 5 keV. In lRussell et al] feOld) . a single-temperature fit to the 
cluster spectrum of Abell 2146 using an absorbed thermal plasma 
emission model yields a temperature of 6.7*2 i keV. The higher X- 
ray temperature measurement is not surprising since M3's derived 
temperature estimate refers to the mean cluster gas temperature 
within r 2 oo and therefore averages over scales where the temper- 
ature is lower. Moreover, emission-weighted temperatures will be 
higher than mass weighted temperature estimates. 



7.3.5 Total mass 

Analytical and numerical simulations have already established 
the integrated SZ signal a s a robust too l for d e termining the 
total c lus ter mass (see e.g. iBartlett & Silk il994h; Barbosa et al.l 
dl996h: lEke, Cole & Frenkl 1 19961) : |Pa Silva et al j fcOOOfc 
iKravtsov, Vikhlinin & Nasal j2006h : iNagail d2006h : iMotl et all 
I2005J)). The measured SZ signal is sensitive to large scales away 
from the cluster core and is therefore able to provide an estimate 
for the Mr which is independent of the small-scale mechanisms 
that regulate the state of the cluster gas near the core. 

We find that, subject to the assumptions of M3 described 
in Section |5~.2.3| at the virial radius, r 2 oo, M T = (4.1 ±0.5) x 
10 14 /r'M Q ; note that this estimate is free from the assumption of 
hydrostatic equilibrium. 



8 CONCLUSION 

The AMI 16-GHz observations of Abell 2146 presented in this pa- 
per show the Sunyaev-Zel'dovich effect produced by this cluster 
with a peak signal-to-noise ratio of 13cr. We detect 15 4-cr„ sources 
within 0.1 of the primary beam in the SA pointed map using the 
high resolution LA observations. These sources were subtracted 
from the SA maps at the LA position using the flux densities ob- 
tained from running our Bayesian analysis software, McAdam, on 
the cluster model M3. Despite the substantial radio emission from 
the direction of Abell 2146, no significant contamination from ra- 
dio sources is visible on the maps. 

We compare our SZ observations with X-ray data taken by 
Chandra and find an offset between the peaks of the two signals. We 
show that the SA data resolves our SZ decrement and note that the 
directions of the most pronounced elongations in the SZ and X-ray 



signals seem to be at =s 90° to each other. These results show com- 
plex dynamics indicative of a cluster merger and the differences in 
the gas emission and pressure distributions. 

We run McAdam, on three different cluster models, all of 
which assume an isothermal, spherical /J-model, and extract pos- 
terior probability distributions of large-scale cluster parameters of 
Abell 2146 in the presence of radio point sources, primordial CMB 
and receiver noise. In Ml, a model representative of more tradi- 
tional cluster parameterizations, the seemingly inconspicuous pri- 
ors on the sampling parameters lead to an effective prior on the de- 
rived parameter f g (r2oo) which biases low this parameter and leads 
to further biases in other model parameters. 

M2 and M3 exploit the observation that the gas fractions do 
not appear to vary greatly between clusters and sample directly 
from /gfooo) — introducing further constraints in our parameter 
space and avoiding the bias problem in Ml . The difference between 
M2 and M3 lies in the derivation of the global cluster gas tempera- 
ture, T. M2 assumes the cluster is in hydrostatic equilibrium while 
in M3 T can be deduced from the virial theorem (assuming all of 
the kinetic energy is in the form of internal gas energy). Given the 
relative masses of the two merging systems in Abell 2146 and, pro- 
vided the primary cluster was virialized before the merger, we find 
that the T derived from the M-T relation in M3 will change by 
» 10% K during the merger. 

The results from M2 and M3 are consistent, despite differ- 
ences in the mean values of the large-scale cluster parameters. 
However, we choose to focus on the results obtained in M3 since 
this model overcomes some of the shortcomings of more tradi- 
tional models and its global temperature estimate is not signifi- 
cantly affected by the merger event. We find that at fioo = 
(4.1 ± 0.5) x 10 14 /i-'M o , = 1.7 ± 0.3, T = 4.5 ± 0.5 keV and 
core radius r c = 358 ± 100/T 1 kpc. We also find that the probability 
of SZ + CMB primordial structure + radio sources + receiver noise 
to CMB + radio sources + receiver noise is 3 x 10 6 : 1. 
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